On the stability of classical chaotic motion under system's perturbations 



o 
o 

(N 
G 



Q 

U 
d 



(N 
> 

m 
O 
O 
oo 
o 

O 



> 

x 



Giuliano Benenti'- 1 *', Giulio CasatiA 1 ' 2 ), and Gregor Veble^ 1 ' 3 ) 
^International Center for the Study of Dynamical Systems, Universita degli Studi dell'Insubria and 
Istituto Nazionale per la Fisica della Materia, Unita di Como, Via Valleggio 11, 22100 Como, Italy 
^ Istituto Nazionale di Fisica Nucleare, Sezione di Milano, Via Celoria 16, 20133 Milano, Italy 
Center for Applied Mathematics and Theoretical Physics, University of Maribor, Krekova ul. 2, SI-2000 Maribor, Slovenia 

(December 5, 2002) 

We study in detail the time behavior of classical fidelity for chaotic systems. We show in particular 
that the asymptotic decay, depending on system dynamical properties, can be either exponential, 
with a rate determined by the gap in the discretized Perron-Frobenius operator, or algebraic, with 
the same power as for correlation functions decay. Therefore the decay of fidelity is strictly connected 
to correlations decay. 



PACS numbers: 05.45.-a, 05.45. Pq 

As it is known, exponential separation of orbits start- 
ing from slightly different initial conditions has been as- 
sociated to classical chaos. It has been remarked that 
the situation in quantum mechanics is drastically dif- 
ferent. Indeed the scalar product of two states (ipi]^) 
is time independent. This has led to the introduction 
of fidelity as a measure of stability of quantum motion 
[1]. More precisely one considers the overlap of two 
states which, starting from the same initial conditions, 
evolve under two slightly different Hamiltonians H and 
H e = H Q + eV. The fidelity is then given by f(t) = 
\(ip\e-xp(iH e t/h)ex_p(-iH t/h)\Tjj}\ 2 . The quantity /(*) 
can be seen as a measure of the so-called Loschmidt echo: 
a state \ip) evolves for a time t under the Hamiltonian Hq, 
then the motion is reversed and evolves back for the same 
time t under the Hamiltonian H e and the overlap with 
the initial state \ip) is considered. 

However, we would like to stress that, in principle, such 
difference between classical and quantum mechanics ac- 
tually does not exist. The Liouville equation, which de- 
scribes classical evolution, is unitary and reversible as 
the Schrodinger equation. However, as stressed in sev- 
eral occasions (see, e.g., Ref. [2]), there exist time scales 
up to which quantum motion can share the properties of 
classical chaotic motion including local exponential in- 
stability. Due to the existence of such time scales, what 
may be different, and indeed it is, is the degree of sta- 
bility of dynamical motion. Indeed, as clearly illustrated 
in the analysis of Loschmidt echo in [3] , quantum motion 
turns out to be more stable than the classical motion. 

The growing interest in quantum computers has at- 
tracted recent interest in this quantity as a measure of 
the stability of quantum computation in the presence of 
hardware imperfections or noisy gates operations. Con- 
fining ourselves to classically chaotic systems, the emerg- 
ing picture which results from analytical and numerical 
investigations [4-12] is that both exponential and Gaus- 
sian decay are present in the time behavior of fidelity. 
The strength of the perturbation determines which of the 



two regimes prevails. The decay rate in the exponential 
regime appears to be dominated either by the classical 
Lyapunov exponent or, according to Fermi golden rule, 
by the spreading width of the local density of states. 

In addition, at least for short times, the decaying be- 
havior depends on the initial state (coherent state, mix- 
ture, etc.). While it can be true that, for practical pur- 
poses, the short time behavior of fidelity may be the most 
interesting one, it is also true, without any doubt, that 
in order to have a clear theoretical understanding and 
identify a possible universal type of quantum decay one 
needs to consider the asymptotic behavior of fidelity. 

This arises the problem of understanding the corre- 
sponding classical decay of fidelity and later on inquiring 
about the times scales at which quantum decay mimics 
the classical one. In the present paper we concentrate 
our attention on the classical behavior. 

What do we know about the decay of classical fidelity 
for chaotic systems? What is the relation with corre- 
lation functions? Can we derive the decay of fidelity 
from the behavior of correlations or is fidelity a com- 
pletely independent function? In a recent paper [12] it 
has been found numerically that, after an initial tran- 
sient, classical fidelity decays exponentially and the rate 
is given by the Lyapunov exponent (see also [8,13]). This 
is in also in agreement with previous papers [4-6] indi- 
cating that quantum fidelity, for strong enough quantum 
perturbation (which, for a fixed classical perturbation 
strength, corresponds to semiclassical region), decays ex- 
ponentially with a rate given by the Lyapunov exponent 
of the corresponding classical system. On the other hand 
we know that the decay of correlation functions is not 
ruled by Lyapunov exponent. In the first place there is 
the general phenomenon of long time tails which means 
power law decay. In addition, for the special cases in 
which one can prove exponential decay, the rate is de- 
termined by the gap in the discretized Perron-Frobenius 
operator and not by the Lyapunov exponent. 

In this paper, we show that the asymptotic decay of 
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classical fidelity for chaotic systems is not related to the 
Lyapunov exponent: Similarly to correlation functions, 
this decay can be either exponential or power law. In the 
first case, the decay rate is determined by the gap in the 
discretized Perron-Frobenius operator, in the latter case 
the power law has the same exponent as for correlation 
functions. 

The classical fidelity f(t) is defined as follows: 

f(t) = [ dxp £ (x,t)po(x,t), (1) 
Jn 

where the integral is extended over the phase space, and 

p (x,i) = t/*p(x,0), p e (x,f) = I^p(x,0) (2) 

give the evolution after t steps of the initial density 
p(x, 0) (assumed to be normalized, i.e. J dxp 2 (x, 0) = 1) 
as determined by the t-th iteration of the Frobenius- 
Perron operators Uq and U e , corresponding to the Hamil- 
tonians Hq and H e , respectively. The above definition 
can be shown to correspond to the classical limit of quan- 
tum fidelity [8,10]. In the ideal case of perfect echo 
(e = 0), the fidelity does not decay, f(t) = 1. How- 
ever, due to chaotic dynamics, when e ^ the classical 
fidelity decay sets in after a time scale 

required to amplify the perturbation up to the size v 
of the initial distribution, with A Lyapunov exponent. 
Thus, for t ^> t v the recovery of initial distribution via 
the imperfect time- reversal procedure fails, and the fi- 
delity decay is determined by the decay of correlations 
for a system which evolves forward in time according to 
the Hamiltonians Ho (up to time t) and H e (from time 
t to time 2t). This is conceptually similar to the "prac- 
tical" irreversibility of chaotic dynamics: due to the ex- 
ponential instability, any amount of numerical error in 
computer simulations rapidly effaces the memory of the 
initial distribution [3]. In the present case, the coarse- 
graining which leads to irreversibility is due not to round- 
off errors but to a perturbation in the Hamiltonian. 

In the following, we illustrate this general phenomenon 
in standard models of classical chaos, characterized 
by uniform exponential instability (the sawtooth map), 
marginal stability (the stadium billard), or mixed phase 
space dynamics (the kicked rotator). 

The sawtooth map is defined by 

p = p + Fo(9), 9 = 9+p, (4) 

where (p, 9) are conjugated action-angle variables, F a — 
K a (9 — 7r), and the overbars denote the variables af- 
ter one map iteration. We consider this map on the 
torus < 9 < 27r, — nL < p < ttL, where L is an 
integer. For K o > the motion is completely chaotic 



and diffusive, with Lyapunov exponent given by A = 
ln{(2 + K Q + [(2 + K ) 2 - 4] 1/2 )/2}. For K > I one 
can estimate the diffusion coefficient D by means of the 
random phase approximation, obtaining D w (n /3)Kq. 
In order to compute the fidelity (1), we choose to perturb 
the kicking strength K = K a + e, with e <C K - In prac- 
tice, we follow the evolution of 10 8 trajectories, which are 
uniformly distributed inside a given phase space region 
of area A n at time t = 0. The fidelity f(t) is given by 
the percentage of trajectories that return back to that re- 
gion after t iterations of the map (4) forward, followed by 
the backward evolution, now with the perturbed strength 
K, in the same time interval t. In order to study the ap- 
proach to equilibrium for fidelity, we consider the quan- 
tity g(t) = (/(*)- /(oo))/(/(0)-/(oo)); in this way g(t) 
drops from 1 to when t goes from to oo. We note that 
/(0) = 1 while, for a chaotic system, /(oo) is given by the 
ratio Aq/A c , with A c the area of the chaotic component 
to which the initial distribution belongs. 
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FIG. 1. Decay of the fidelity g(t) for the sawtooth map 
with the parameters K — (\/5 + l)/2 and e = 1CP 3 for dif- 
ferent values of L = 1, 3, 5, 7, 10, 20, oo from the fastest to the 
slowest decaying curve, respectively. The initial phase space 
density is chosen as the characteristic function on the support 
given by the (6,p) € [0, 2ir) x [-tt/100, tt/100]. Note that be- 
tween the Lyapunov decay and the exponential asymptotic 
decay there is a oc l/Vt decay, as expected from diffusive be- 
havior. Inset: magnification of the same plot for short times, 
with the corresponding Lyapunov decay indicated as a thick 
dashed line. 

The behavior of g(t) is shown in Fig. 1, for Ko = 
(V5 + l)/2 and different L values. One can see that only 
the short time decay is determined by the Lyapunov ex- 
ponent. It takes place for t v < t < t e , with t v defined 
in Eq. (3) and t e = (1/A) ln(27r/e) time scale required to 
amplify the effect of the Hamiltonian perturbation up to 
the maximum extension in the angle 9. The Lyapunov 
regime is followed by a power law decay [12] oc 1/VDt 
until the diffusion time to ~ L 2 /D and then the asymp- 
totic relaxation to equilibrium takes place exponentially, 
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with a decay rate 7 (shown in Fig. 2), which, as discussed 
below, is ruled not by the Lyapunov exponent but by the 
largest Ruelle-Pollicott resonance [14]. In particular, it 
is e- independent. 




FIG. 2. Asymptotic exponential decay rates of fidelity for 
the sawtooth map (K = (V5 + l)/2, e = 1CT 3 ) as a function 
of L. Fhe rates are extracted by fitting the tails of the fidelity 
decay in the Fig. 1 (triangles) and from the discretized Per- 
ron- Frobenius operator (circles) . The line denotes the oc 1 / L 
behavior of the decay rates, predicted by the Fokker-Planck 
equation, which describes the classical motion in the diffusive 
regime. 



We determine numerically these resonances for the 
sawtooth map using the following method [15,16]: 

(i) the phase space torus (0 < 9 < 2tt, —irL < p < ttL) 
is divided in N x NL square cells; 

(ii) the transition matrix elements between cells are de- 
termined numerically by iterating for one map step the 
phase space distributions given by the characteristic func- 
tions of each cell: in this way we build a finite dimensional 
approximation of the one-period evolution operator Uo; 

(iii) this truncated evolution matrix Uq N ^ (of size LN 2 x 
LN 2 ) is diagonalized: it is no longer unitary, and its 
eigenvalues z[ N ^ are inside the unit circle in the complex 
plane (an example is shown in Fig. 3). The non-unitarity 
of the coarse-grained evolution is due to the fact that 
the transfer of probability to finer scale structures in the 
phase space is cut-off, and this results in an effective dis- 
sipation [15]; 

(iv) resonances correspond to "frozen" non-unimodular 
eigenvalues, namely z\ N ^ — > when N — > 00, with 
\zi\ < 1. Convergence of eigenvalues to values inside the 
unit circles comes from the asymptotic self-similarity of 
chaotic dynamics [15]. 

As it is known, the asymptotic (t — > 00) relaxation of 
correlations is determined by the resonance with largest 
modulus, \z\ = maxi|zi| < 1, giving a decay rate 70 = 
ln\z\. In Fig. 2 we illustrate the good agreement between 
the asymptotic decay rate of fidelity (extracted from the 



data of Fig. 1) and the decay rate as predicted by the gap 
in the discretized Perron-Frobenius spectrum. It should 
be stressed that, since fidelity involves forward and back- 
ward evolution, the fidelity decay at time t has to be 
compared with the correlations decay at time 2t. For 
this reason in Fig. 2 the circles correspond to 7 = 270. 
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FIG. 3. Spectrum of the discretized Perron-Frobenius oper- 
ator for the sawtooth map with parameters Ko = (\/5 + l)/2, 
L = 7 and discretization N = 20. The asymptotic decay of 
fidelity is determined by the largest modulus eigenvalue apart 
from the eigenvalue 1. 
^0. 




FIG. 4. Power law decay of fidelity for the stadium billiard 
with radius R — 1 and length of the straight segments do = 2 
(the perturbed stadium has d = do+e, with e = 2x 1CF 3 ). The 
initial phase space density was chosen to be a direct product 
of a characteristic function on a circle in configuration space, 
the center of which was at (0.5,0.25) as measured from the 
center of the billiard and its radius was 0.1, while for mo- 
menta the <5(|p| — 1) distribution was used. The dashed line 
has slope — 1. 



We would like to stress that the same qualitative be- 
havior of Fig. 1 is obtained in the presence of stochastic 
noise, e.g. the backward evolution is driven by a kicking 
strength K(t) = K + e(t), with {e(t)} t =i,2,... uniformly 
and randomly distributed inside the interval [— e, e]. In 
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particular, we observed the initial Lyapunov decay and 
the asymptotic exponential relaxation with the same rate 
7. This means that the effect of a noisy environment on 
the decay of fidelity for a classical chaotic system is sim- 
ilar to that of a generic static Hamiltonian perturbation. 




t 

FIG. 5. Decay of fidelity for the kicked rotator with 
K = 2.5, L = 1, and e = 1CT 3 (full curve). The support of 
the initial (characteristic) density is inside the chaotic compo- 
nent, with (0,p) G [0,0.2] x [0,0.2]. The dotted curve repre- 
sents the exponential decay at a rate given by the Lyapunov 
exponent A « 0.534. The dashed line has slope —0.55. The 
dot-dashed curve gives the correlation decay D(2t), for the 
same initial density and for twice the time t. It is clearly seen 
that, asymptotically, fidelity and correlations have the same 
power law decay oc t~ 0,55 . 

Further confirmation for the validity of the above illus- 
trated scenario has been obtained by analyzing systems 
in which the asymptotic decay of correlations is algebraic. 
This happens in two cases: 

(i) when the system possesses marginally stable orbits. 
A typical example is the stadium billiard in which, as it 
is known [17], correlations decay as 1/t] 

(ii) when there is mixed phase space [18]. A typical ex- 
ample is the kicked rotator model (described by Eq. (4) 
with F — Kq sin 9) . 

Since in the long time limit the fidelity decay at time t is 
still related to the decay of correlations at time 2t, in the 
case of power law decay of correlations as t~ a , we expect 
a power law decay of fidelity with the same exponent a. 
This is indeed confirmed by our numerical results. In 
Fig. 4 it is shown that, for the stadium billiard, fidelity 
decays asymptotically oc 1/t, as expected. In Fig. 5 
we compare the fidelity decay (at time t) and the cor- 
relations decay (at time 2t) for the kicked rotator with 
kicking parameter K — 2.5, for which the phase space 
contains chaotic components and stable islands. The cor- 
relator is given by D(t) = (C(t)-C(oo))/(C(0)-C(oo)), 
with C(t) = J n dxpo(x,t)p(x,0). It is seen that, after 
an initial Lyapunov decay, fidelity approaches the same 



asymptotic power law decay of correlations [19]. 

In summary, we have shown that in chaotic systems the 
asymptotic decay of classical fidelity, which describes the 
structural stability of motion under system's perturba- 
tions, is analogous to the asymptotic decay of correlation 
functions. This asymptotic decay can be either exponen- 
tial or algebraic, depending on the dynamical properties 
of the system. In any instance, it is not related to the 
local exponential instability ruled by the Lyapunov ex- 
ponent and it is e- independent. It would be interesting 
to understand what are the implications of these findings 
for the decay of quantum fidelity. 
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